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ABSTRACT 

Recent quasar microlensing observations have constrained the X-ray emis- 
sion sizes of quasars to be about 10 gravitational radii, one order of magnitude 
smaller than the optical emission sizes. Using a new ray-tracing code for the 
Kerr space-time, we find that the observed X-ray flux is strongly influenced by 
the gravity field of the central black hole, even for observers at moderate incli- 
nation angles. We calculate inclination-dependent flux profiles of active galactic 
nuclei in the optical and X-ray bands by combining the Kerr lensing and pro- 
jection effects for future references. We further study the dependence of the 
X-ray-to-optical flux ratio on the inclination angle caused by differential lensing 
distortion of the X-ray and optical emission, assuming several corona geometries. 
The strong lensing X-ray-to-optical magnification ratio can change by a factor 
of ~10 for normal quasars in some cases, and another factor of ~10 for broad 
absorption line quasars (BALs) and obscured quasars. Comparing our results 
with the observed distributions in normal and broad absorption line quasars, we 
find that the inclination angle dependence of the magnification ratios can change 
the X-ray-to-optical flux ratio distributions significantly. In particular, the mean 
value of the spectrum slope parameter a ox , 0.3838 log F 2 k e v/-^25ooA! can differ 
by ~0. 1-0.2 between normal and broad absorption line quasars, depending on 
corona geometries, suggesting larger intrinsic absorptions in BALs. 

Subject headings: accretion, accretion disks — black hole physics — gravitational 
lensing — quasars 
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1. Introduction 

X-ray emission is a defining characteristic of AGN; however, unlike optical emission, 
its origin is currently unclear. In contrast to X-ray binaries, standard AGN accretion disk 
theory predicts disk temperatures too low to emit X-rays. The X-ray emission is thought 
by many to be generated by inverse Compton scattering of UV photons emitted from the 
disk by hot electrons in a corona above/around the disk (see review of Reynolds & Nowak 
2003). However, both the geometry (e.g., its size, shape, and position, e.g., height above 
the accretion disk) and the physics (e.g., the mechanism for heating the electrons) of the 
corona are under debate (Haardt & Maraschi 1991; Haardt et al. 1994; Collin et al. 2003; 
Torricelli-Ciamponi et al. 2005; Trzesniewski et al. 2011). Based on the Bayesian Monte Carlo 
microlensing analysis technique (Kochanek 2004; Kochanek et al. 2007), recent observations 
in gravitational microlensing have constrained the X-ray emission size of quasars to be of 
order 10 r g (gravitational radius, GM/c 2 , Blackburne et al. 2006; Pooley et al. 2006; Morgan 
et al. 2008; Chartas et al. 2009; Dai et al. 2010; Blackburne et al. 2011), about one order 
of magnitude smaller than the optical emission region (Poindexter et al. 2008; Morgan et 
al. 2010; Mosquera et al. 2011). Chen et al. (2011, 2012a) detected energy-dependent X-ray 
microlensing in Q 2237+0305 and several other lenses for the first time, and their results 
suggest that the hard X-ray emission might come from regions smaller than the soft, while 
in RXJ 1131-1231 the effect is the opposite (Chartas et al. 2012). Morgan et al. (2012) 
found that the majority of the X-ray emission in QJ 0158—4325 comes from regions close to 
6r g , i.e., in the vicinity of the last stable orbit of a Schwarzschild black hole. 

For distances less than a few gravitational radii, the flux and profile of both the X-ray 
continuum and metal emission lines (in particular, the FeKct line) are strongly influenced 
by the gravitational field of the supermassive black hole at the center of an AGN. Due to 
the complexity of the Kerr metric (space-time of a black hole with angular momentum, or 
spin), much work has focused on either the analytic treatment of the geodesies in a Kerr 
metric (Kraniotis 2004, 2005; Hackmann et al. 2010) or on developing fast ray-tracing codes 
(Karas et al. 1992; Dexter & Agol 2009; Vincent et al. 2011). There has also been much 
work studying AGN X-ray spectra under strong gravity (Fabian et al. 1989; Chen et al. 
1989; Laor et al. 1989; Laor 1991; Hameury et al. 1994; Bromley et al. 1997; Fabian & 
Vaughan 2003; Beckwith & Done 2004; Fuerst & Wu 2004; Schnittman & Krolik 2010). 
The observed broad, skewed iron lines (e.g., MCG6— 30— 15, Tanaka et al. 1995) strongly 
suggest that the reflection component near the black hole contributes significantly to the 
total emissivity. Since the iron line profile strongly depends on the black hole spin, and 
the relativistic distortion of the iron emission line is much more easily measured than the 
broad X-ray continuum, most of these papers concentrated on the iron emission line. For 
example, Fuerst & Wu (2004) studied the radiation transfer of emission lines in the Kerr 
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space-time. However, less work has been done on the continuum emission. It has been found 
that strong lensing of an accretion disk by the central supermassive black hole is strongly 
dependent on inclination angles (Cunningham 1975), and for very large inclination angles 
higher order images can contribute a significant fraction of the total luminosity (Beckwith 
& Done 2005). Polarization of the continuum emission near black holes was investigated in 
Connors et al. (1980) and is becoming an important tool in measuring the inclination angle 
of an accretion disk and the spin of its black hole (Dovciak et al. 2004; McClintock et al. 
2006; Li et al. 2009; Schnittman & Krolik 2009; Dovciak et al. 2008, 2011). Quasi-periodic 
oscillations (QPOs) have been observed in the X-ray light curves of black hole binaries and 
black hole candidates for more than a decade (Remillard & McClintock 2006); however, no 
convincing QPO signal has been detected for AGN until very recently (Gierlinski et al. 2008; 
Middleton et al. 2009). High frequency QPOs are often modeled in the context of strong 
gravity (the oscillation period is comparable to the Keplerian orbital period near the inner 
most stable orbit), and this can be used to accurately measure the black hole spin once 
the correct model of QPOs has been found (Wagoner et al. 2001; Schnittman 2005). The 
study of radiation transfer is crucial for these subjects, and Monte Carlo ray-tracing methods 
have been used extensively to solve (continuum or line) radiation transfer problems in AGNs 
(Agol & Blaes 1996; Laurent & Titarchuk 1999; Bottcher et al. 2003; Dolence et al. 2009; 
Schnittman & Krolik 2010). However, modeling of the AGN X-ray continuum emission was 
hampered by the lack of an accurate estimate of the size of the emission region. We study 
the influence of the strong field of the massive spinning black hole centering the AGN on 
the X-ray continuum emission using the most recent constraints from quasar microlensing 
observations and our ray-tracing code. Despite the fact that there are nicely written public 
codes for ray-tracing in Kerr space-time available (e.g., Dexter & Agol 2009), we develop our 
own ray-tracing code (specially designed for backward ray-tracing, see §3.2[) with a future 
goal to solve the radiation transfer equation in Kerr space-time using the characteristics 
method and operator perturbation algorithm (Hauschildt & Baron 2006; Chen et al. 2007; 
Baron et al. 2009a,b). 

In this paper, we study the inclination angle dependence of the flux profiles in the X-ray 
and the optical continuum of AGNs and their ratios caused by strong gravity, providing a 
reference for future usage. The X-ray-to-optical flux ratio, conventionally characterized by 
a parameter a ox (Tananbaum et al. 1979), is empirically measured in many AGNs. Our 
inclination dependent a ox with the Kerr lensing and projection effects taken into account 
will provide a baseline model for interpreting the observed a ox in AGNs. For example, 
assuming that quasars are observed at randomly distributed solid angles, if there exists 
strong inclination angle dependence for the X-ray-to-optical flux ratios caused by the strong 
gravity field, the flux ratios will be observed with larger fluctuations than the intrinsic ones. 
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The mean value, scatter, and shape of the distributions might be changed by the strong 
gravity (see $6]) . Therefore, this can provide a test for different corona geometries, combined 
with a large number of measurements of a ox from different AGNs. 

Based on the unification model (Antonucci 1993), AGN can be classified as Type I 
(normal) or Type II (obscured, possibly by a torus) with ~50% Type II AGN locally at 
low redshifts (e.g., Simpson 2005; Maiolino et al. 2007; Gilli et al. 2007). There are sugges- 
tions that the fraction of Type II AGN decreases with luminosity or increases with redshift 
(Hasinger 2008; Treister et al. 2009; Bianchi et al. 2012). Because emission of Type II AGN 
is strongly modified by Compton-thick absorption, which makes it difficult to study the in- 
trinsic a ox , we do not compare our model predictions to Type II AGN. In Type I quasars, a 
significant fraction (20-40%) are broad absorption line (BAL) quasars (e.g., Dai et al. 2008; 
Shankar et al. 2008). These BALs are only moderately obscured, and several studies show 
that their intrinsic a ox are consistent with non-BALs (Green & Mathur 1996; Brandt et al. 
2000; Green et al. 2001; Gallagher et al. 2002). Therefore, we compare our model predictions 
with a ox measurements in Type I quasars including BALs. If the X-ray-to-optical magnifi- 
cation ratio Rx/o is different between normal quasars and BALs, then there might be bias 
in the current estimate of BAL absorption caused by ignoring the effect of strong gravity. 
We investigate these ideas in this paper. 

In §2] we introduce the ray-tracing code for Kerr space-time developed in this paper. 
In §3] we present the basic formalism for studying strong lensing of Kerr space-time. In §H 
we study the influence of strong gravity on the optical emission. In §0 we study strong 
lensing of the X-ray emission assuming three X-ray corona geometries. In £0we compare the 
lensed a ox of normal/BAL quasars with observations, and estimate possible errors in current 
measurements of BAL absorption caused by ignoring effects of the Kerr strong gravity. In 
§[7] we draw our conclusions. Readers not interested in mathematical details are invited to 
skip §0 §E21 and all or part of g5J 



2. Ray- Tracing In Kerr Space-time 

A Kerr black hole is fully characterized by two parameters, i.e., mass M and angular 
momentum a. In Boyer-Lindquist coordinates (Boyer & Lindquist 1967), the Kerr metric 
(Kerr 1963) with signature (— , +, +, +) can be written as 

v2 2 

ds 2 = -a 2 dt 2 + — sin 2 9(d<P - uodtf + ^-dr 2 + p 2 d9 2 (1) 
p 2 A 
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where t, r, 9, and (p are the four independent coordinate variables, and 



P = 

A = 

£ 2 = 

a 2 = 

W EE 
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(2) 



The off-diagonal terms of the Kerr metric are characterized by u, which represents the 
dragging of the inertial frame. A photon's geodesic motion can be described by a Hamiltonian 
T-L(x a ,p a ) defined on an eight-dimensional phase space (the index a runs from to 3), 
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2p 2 Asin 2 # ^' 
where 

g ab ig 

the inverse metric tensor and p a is the covariant photon 4-momentum vector 
in the coordinate basis. To obtain photon trajectories in the Kerr space-time, we solve the 
group of eight first order Hamilton equations 



dx a 

dp g 
di 



dU 

dp a 

&H 

dx a 



(4) 



where ^ is the affine parameter (refer to Chandrasekhar 1983 for details), instead of the four 
second order geodesic equations 



d 2 x a dx b dx c 

+ F 6c— — = °> 



(5) 



de oc di di 

where T a bc are the Christoffel symbols. The number of equations reduces to six since p t and 
are motion constants, and we give them in the Appendix. These six first order differen- 
tial equations are much easier to solve numerically than the four second order differential 
equations. There is another motion constant discovered by Carter (1968), 

,,2 
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+ cos 2 9 



sin 2 6 



2 2 

a p t 



(6) 



which we use as an independent check for our ray-tracing code (the maximum relative error 
of C allowed in the ray-tracing is less than 10~ 7 ). We solve Eqs. fl4]) numerically using a fifth- 
order Runge-Kutta algorithm with adaptive step size control (see e.g., Press et al. 1992). 
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Let u and u e be the photon frequency measured by a distance observer and an observer 
co-moving with the light source. We define the redshift factor 

_ v Q (-u a p a ) 
9 = — = 7 — ;r~T' ( 7 ) 

U e {-U a p a ) e 

where and «" are the four-velocity of the distant observer and the source. For the distant 
observer, we choose 

u a = -(1,0,0, u), u a = (-a, 0, 0, 0) (8) 
a 

(the so-called zero-angular-momentum observer). An advantage of this choice is that the 
collection {u a , d r , do, d^} forms an orthogonal basis and therefore the photon's spatial di- 
rection can be easily specified with respect to {d r , do, d^}, whereas <9 r , <9^} does not 
form an orthogonal basis (Bardeen et al. 1972). The four- velocity of the light emitter is 
model-dependent. We assume that the source is circling around the z-axis, i.e., 

M a = 7 (l,0,0,fi). (9) 

The selection of Q is not arbitrary, but constrained by the causality condition u a u a = —1. If 
the source is on the accretion disk (e.g., optical emission), we assume a prograde^ Keplerian 
flow for the fluid (Bardeen et al. 1972), i.e., 

" = n* - r 3/* + aM v3 - < 10 > 

If the source is not in the equatorial plane (e.g., the X-ray ball model in §5.21) . we use the 
profile of Ruszkowski & Fabian (2000), 



9 

V/2 



u, (11) 



where 9 is the spherical polar angle in the Boyer-Lindquist coordinates (an interpolation 
between fl K and uj of the frame dragging, see Eq. (j2J)). 



3. Basics of Strong Lensing in Kerr Space-Time 

Strong lensing effects of a Kerr black hole on the accretion disk and X-ray corona differ 
from standard linear gravitational lensing theory in a few important aspects. First, linear 



1 Here prograde means that a • J > where a is the (vector) spin of the black hole, and J is the angular 
momentum of the accretion disk. 
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lensing theory requires a very small bending angle, i.e., GM/tq <C 1 where M and ro are the 
lens mass and the impact vector, respectively. This assumption breaks down for accretion 
disks whose inner cutoff can be as close as a few gravitational radii from the lensing black 
hole. Second, the source, lens, and observer are far away from each other in the linear lensing 
theory, whereas for Kerr lensing, there is no well-separated lens and source planes, i.e., the 
source is very close to the lens. Third, the space-time in the source plane is assumed to be flat 
in the linear theory, whereas for Kerr lensing, this is not true. Last, since both the source and 
the observer are far away from the lens in the linear theory, there is no gravitational redshift 
caused by the lens, and the redshift comes solely from cosmology. This is incorrect for Kerr 
lensing. Since photons are emitted by sources moving with relativistic speeds in strongly 
curved space-time, both the Doppler and gravitational redshifts have to be taken into account 
in Kerr lensing. Consequently, there is no simple lens equation relating the source and image 
positions through a single bending angle (the bending happens continuously, instead of at 
a single point in the lens plane as assumed in the linear lensing theory). We identify the 
source position for each ray arriving at the observer through backward ray-tracing in the 
Kerr space-time (see §3.2|) . The redshift effect is very important for computing observable 
quantities from the accretion disk around the supermassive black hole. 

By definition, the monochromatic flux F Va measured by an observer O sitting in a 
(asymptotically) flat space-time is 



where dVL Q is the differential solid angle measured at the observer, \i is the angle cosine with 
respect to the normal of the detector window (see Figured]), and can be safely dropped for 
sources very far away from the observer, i.e., the source is observed with a tiny solid angle. 
In standard lensing theory, there are no frequency redshifts between the source and observer 
(besides the cosmological redshift), therefore I„ is a conserved quantity along the light path 
in case of no scattering and absorption, thus 



where 6 and (3 are the image and source positions, respectively. Since the space-time in 
both the source and lens plane are assumed to be flat, the integral over solid angle can be 
easily transformed into 2-D surface integrals in the lens or source plane. To compute the 
flux magnification, we simply compute the ratio of the image area to the source area in the 
lens plane, i.e., 






(13) 




(14) 
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This is equivalent to computing the ratios of the solid angles subtended by the image and 
the source at the observer. On the other hand, for Kerr lensing, we have 

F Vo = J I° o bs (n)d0 o = J g 3 (n)I^ cc (x a (n),p a (h))dQ , (15) 

where n is the photon's 3-D direction at the observer, (x a ,p a ) is the photon's position and 
momentum at the emitter (found through backward ray-tracing), u e is the source frequency, 
and g is the redshift factor defined earlier. Clearly, the effect of Kerr lensing strongly depends 
on the redshift factor g, and consequently on the four- velocity of the source. Solid angles 
(at the observer) and areas (in the lens or source plane) are often used interchangeably in 
linear lensing theory; however, in Kerr lensing theory, solid angles at the observer are more 
appropriate quantities to use since the space-time around the source and lens is strongly 
curved E 



3.1. Models of Optical and X-Ray Emission 



We assume that the optical/UV emission comes from the accretion disk (Shakura & 
Sunyaev 1973), and is Planckian, 



B(u,r) 



2h 



cxp 



hu 
kT(r) 



(16) 



The effective temperature can be approximated as T(r) oc r -3 / 4 ^ 1 / 4 where k = 1 — a/ ?"isco / r 
(risco is the innermost stable circular orbit, Bardeen et al. 1972). For convenience we set 
k to unity since this factor only influences the innermost region of the accretion disk which 
does not contribute significantly to the total optical flux. Likewise we have ignored the 
general relativistic corrections to T(r) in the innermost regions (Novikov & Thorne 1973). 
From Eq. ffl5|) . the observed monochromatic flux F Vo at frequency v Q is 



g 3 B(v e ,r)dtt 



2h 



2h 
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9T(r) 
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,4/3 



)dQ a . 



(17) 



2 The familiar angular diameter distance defined in the FLRW cosmology, e.g.,Z? s of the source, can not 
be used to compute the source area directly since it does not respect the space-time curvature at the source 
caused by the supermassive black hole. 
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In the above derivations we have defined the so-called scale radius r s (not the Schwarzschild 
radius) by hv = kT Q = kT{r s ). We assume that the accretion disk moves with the Keplerian 
flow, see Eq. f lTU]) . We choose the inner cutoff of the optical emission disk at the innermost 
stable circular orbit, r ISC o, where r ISCO = 6r 9 , 4.23 r g , and 1.24 r g for a = 0, 0.5M, and 
0.998M black holes, respectively (Bardeen et al. 1972). This last value corresponds to an 
extreme Kerr black hole (Thorne 1974). 

The spectrum of the X-ray emission of AGN is observed to follow a power-law. We 
therefore assume a simple power-law model for the X-ray emission, 

^MoclJ^, (18) 

where r is the radial coordinate of the source, n specifies the steepness of the radial profile, 
and T is the photon index. Using Eq. ( I18p . we also assume that the X-ray emission is isotropic 
in the local rest frame of the X- rayso^eS The observed monochromatic flux is 

F ^ = 1 9 Z hJ^o = [ p^zrda, = / ^T^o = ^ I (19) 

J J e J o o J 

which follows a power-law with the same photon index V and we have dropped the unim- 
portant constant in Eq. ( Tl8l) . The strong lensing correction is contained in the last integral 
in the equation above, and the redshift factor g depends only on the photon trajectory, or 
equivalently on the photon's arriving direction n at the observer, but not on its frequency 
because gravitational lensing is achromatic. Different radial profiles and sizes of X-ray emis- 
sion regions give different strong lensing corrections. The observed bolometric flux (e.g., 
1-50 keV, in the observer's frame) is 

/* l^m ax p fm ax 1 

F= F Vo dv = / —dv x / y —<m o . (20) 

Since the first integral in the above equation is independent of the second one, the strong 
lensing corrections to the integrated flux and the monochromatic flux are the same for the 
simple X-ray emission profiles considered in this paper. 

We consider three different X-ray geometries: a thin X-ray disk immediately above the 
accretion disk moving with Keplerian flow, see Figure |2] (low X-ray disk model), a spherical 
corona around the black hole (X-ray ball model), and a thin disk above the accretion disk 
with different heights (high X-ray disk model). For the low X-ray disk model, we choose 



3 The isotropy assumption is made for simplicity. However, the formalism developed in this paper does 
not require isotropic intensity profiles. 
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the inner cutoff at risco- We choose no inner cutoff for the high X-ray disk model. For the 
X-ray ball and high X-ray disk model, we assume that the X-ray source is moving with the 
Ruszkowski & Fabian (2000) profile, Eq. f lTTj) . see Figures [3] and HI We take the photon index 
T = 2, which is normal for quasar X-ray spectra. For each X-ray geometry, we evaluate three 
spin parameters: a = 0, 0.5 M, and 0.998 M. 



3.2. Backward Ray- Tracing 

For flat space-time, light paths are straight lines, and there is no gravitational redshift. 
Ignoring the special relativistic Doppler shifts and aberration effect, we have g — 1. Therefore, 
Eqs (fTTI) and (fT9j) can be easily carried out as a surface integral over the accretion disk or 
the X-ray source, i.e., 

F 0( flat ) = 27rCOS0 I B{Vo r y dr ^ (21) 
r obs 



and 



px(flat) 1 / 1 1 f dAj: 

v° -,r-i / r n „r-i~,2 / _ n ' 



K .1 r" K r ohs 



where 9 is the inclination angle of the disk, r Q b s is the radial coordinate of the observer, 
and dAj; is the differential (X-ray) source area orthogonal to the line of sight. For curved 
space-time with relativistic flows, Eqs. ( TT7]) and f lT9~j) need to be computed numerically. 

We choose a pencil beam with a rectangular cross-section, large enough to contain 
all rays arriving at the observer (Figured]), around the central light ray toward the black 
hole, and divide the solid angle space into a uniform grid. Each pixel, labeled by 
corresponds to an image area dA it j = dO^r^,. Since we have used a uniform rectangular 
grid, dQij = dQ is constant for each pixel. Therefore, to compute the image area, we need 
only count the number of pixels whose light rays end at the accretion disk. To compute 
the observed flux, we weight pixels by the integrands in Eqs (fT7]) and f lT9~j) for optical and 
X-ray emission, respectively. Eqs f lT7|) and f[T9~l) are written as integrals over solid angle at 
the observer; however, we emphasize that the integrands are evaluated at the emitter (near 
the black hole), and each piece of solid angle dQ (h) is mapped to a photon's initial phase 
space coordinate (x a ,p a ) on the accretion disk through backward ray-tracing. For example, 
if a photon from the optical accretion disk arrives at the observer along direction n, we find 
its r coordinate on the accretion disk and redshift factor g through backward ray-tracing 
(Figured]). If g < 1 (redshifted), then we sample the temperature of the blackbody radiation 
at radius r/g 4 ^ 3 (> r, T(r/g 4 ^ 3 ) < T(r)). This approach takes into account gravitational 
light deflection, frequency redshift, and area distortion automatically. We do not use the 
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popular transfer function method (Cunningham 1975), which only considers light emitted 
on the accretion disk (the equatorial plane), since we allow X-rays to be emitted outside of 
the equatorial plane. 



3.3. Strong Lensing Correction to the X-Ray To Optical Flux Ratio 

Let Mq and Mx be the strong lensing magnifications of the optical and X-ray emission, 
then the ratio between the observed monochromatic flux in the optical and X-ray band can 
be written 

pX » r T7.X(flat) pX 
t v a _ M X £v _ p(lens) R (proj) f\ 

F9 Mo P°( flat ) ~ x /° x /° F°, 
where we have defined the lensing magnification ratio, 



(fiat) 



^x/o -po 



(flat) 



(23) 



„(lens) _ M X ( . 

Rx/o = Wq > (24) 

the (Euclidean) projection area ratio Rx/o ( W1 th respect to the face on case), and the total 
ratio 

p — p( lens ) p(P r °j) (nr\ 

Kx/o = tt x /o K x/o y Zb > 

which contains all inclination dependence (lensing plus projection). If the X-ray and optical 
emission regions are different in geometrical size, surface brightness profile, or are located 
at different places with respect to the central black hole, strong gravitational lensing of 
the Kerr black hole should magnify or demagnify these two emission sources differently. 
In particular, the inclination angle dependence of the magnification should be different for 
X-ray and optical emission, and the magnification ratio Rx/o should be a function of 6. 
Among the three X-ray geometries considered, Rx/o = 1 101 the low and high X-ray disk 
model because the factor cos# is canceled when considering both the optical and X-ray 
emission, and Rx/o purely reflects the differential lensing effect. For the X-ray ball model, 
Rjuq — (1 + cos6')/2cos^ is inclination-dependent (see §5. 2[) . where the projected X-ray 
emission area is non-vanishing at large inclination angles. Thus, Rx/o is affected by both 
lensing and projection effects for the X-ray ball model. 



We investigate strong lensing of optical and X-ray emission by Kerr black holes in the 

(lens) 
X/O 



next two sections, and pay particular attention to the lensing magnification ratios R^qIO). 



In §|6]we compare the strong lensed modulated a ox with observations. 
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4. Strong Lensing of the Optical Emission 

Because of uncertainties in parameters such as the black hole mass and the accretion ef- 
ficiency, the scale radius r s , or the half light radius (i.e., radius of the circular region centering 
the peak brightness of the emission and containing half the total flux at a given wavelength) 
r haif = 2.45 r s , cannot be precisely determined from the accretion disk theory. Here, we use 
the empirical results from recent quasar microlensing observations (Kochanek 2004; Pooley 
et al. 2006; Morgan et al. 2008; Eigenbrod et al. 2008; Chartas et al. 2009, Dai et al. 2010, 
Blackburne et al. 2011; Blackburne et al. 2012), and assume that rhaif = 100 r g , and take 
r s = 41 r g . We choose the outer boundary of the optical disk at r disk = 200 r g which is 5 
times the scale radius r s , and contains ~80% of the monochromatic flux. We ignore the 
lensing correction for regions with r > 200 r g , and assume that the flux magnification ratio 
is 1 beyond rdisk- Consequently, if the flux magnification is M <2 oor g within rdisk, the total 
flux magnification Mo is then Mo = 0.2 + 0.8M <2 oor g - 

We plot the strong lensing distortion to the image area and specific flux F Vo in Figure 
for three different spin parameters, a = 0, 0.5M, and 0.998M. Since the size of the optical 
emission region is relatively large (of an order 100 r g ), and the intensity profile (blackbody 
radiation with T oc r~ 3//4 , see Eq. (flBl) ). is not very steep, the distortion in either the image 
area or flux is not significant (close to 1) for AGNs observed nearly face on, and the area/flux 
magnification curves Mq{9) are flat until the inclination angle 9 becomes large (> 75°, i.e., 
observed near edge on). Since the black hole spin a contributes only in higher orders (e.g., 
u ~ oM/r 3 ), we do not see a significant difference in the flux or image area distortion 
for different spin parameters for optical emission. If we increase the disk size rd; sk , the 
magnification of image area Mq^ at all inclination angles will eventually approach one, but 
the magnification in flux Mq ux changes only slightly, since the majority of the emission is 
enclosed within rdisk and the effect of gravity is weak outside of rdi sk . 



5. Strong Lensing of the X-Ray Emission 

In this section we consider three different X-ray models: a thin X-ray disk immediately 
above the accretion disk (the low X-ray disk model), a spherical corona around the black hole 
(the X-ray ball model), and a thin disk above the accretion disk with different heights (the 
high X-ray disk model). The results for these X-ray models as well as the optical emission 
model discussed in §4] are summarized in Table [TJ 
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5.1. Low X-Ray Disk Model 

We first consider a simple model where X-rays are emitted by a tenuous disk immediately 
above the quasar accretion disk and this X-ray disk moves with the same Keplerian flow as 
the optical emission disk. This is partially motivated by the so-called "sandwich" corona 
model (Haardt & Maraschi 1991, 1993; Ruszkowski & Fabian 2000). We choose the outer 
and inner cutoff at rdisk = 10r 9 and r in = risco, respectively. For the intensity profile, we 
use Eq. (II 8p . with n = 3 (i.e., l v oc r -3 , following the optical profile) or n = (no radial 
dependence). For the case of an extreme Kerr black hole (a = 0.998 M), we also tried the 
case where r^sk = 5r 9 , which is not realistic for a = or 0.5 M, considering the inner cutoff 
r isco we choose. We show the lensing distorted and redshifted images of such an X-ray disk 
in Figure |6] for observers at 9 — 15°, 45°, and 75° (from nearly face on to nearly edge on), 
and three spins, a = 0,0.5, and 0.998 M. All images in Figure |6] are distorted as results of 
the lensing effect. Disks viewed at large inclination angles show stronger image distortion, 
and black holes with larger spins produce larger redshifts or blue-shifts. Since the optical 
accretion disk (§UJ) is moving with the same Keplerian flow, the images for an optical disk 
are very similar, except that the image sizes are larger. 

In Figure [7] we show the amplification of the X-ray image area and its ratio with respect 
to the optical image area for each case. The curves in the two panels resemble each other 
at 9 < 75° since the optical area magnifications, Mq(9), are nearly flat and very close to 
1 in this range of inclination angles (Figure |5]). At larger inclination angles, the optical 
disk area is also significantly magnified (Figure |5]), and thus the corresponding curves in the 
two panels deviate from each other. The X-ray area amplification increases with inclination 
angles because more photons pass near the black hole at larger inclination angles, which 
are more strongly affected by lensing. The area amplification ratios also increases with 
inclination angles because the X-ray emission region is much smaller, resulting in stronger 
lensing. 

In Figure [S] we show the magnification of the observed flux F v as a function of in- 
clination angle. For small inclination angles, the Doppler effect is unimportant, and the 
gravitational redshift reduces the observed specific flux through the invariance of I v jift . For 
large inclination angles, the amplification of image area increases rapidly (with 9) because 
of strong lensing (Figure [7]), and Doppler blue-shifted regions (approaching the observer) 
compete with Doppler redshifted regions (receding from the observer). The lensed specific 
flux begins to be greater than the unlensed flux, i.e.,Mx > 1 at 9 > 70°. The model clearly 
predicts an inclination angle dependence of the magnification ratio Rx/o- Furthermore, the 
dependence of Rx/o on inclination angles is stronger for smaller X-ray disks, larger spins 
(smaller inner cutoff risco)> or steeper radial profiles (more weight on the central part close 
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to the black hole). If the sizes of the X-ray and optical emission regions constrained from 
microlensing observations are correct, then even the weakest inclination angle dependence 
among our models, i.e., case in which a = 0, r^ s ^ = 10 r g , and flat radial profile, gives an 
increase of a factor of two in Rx/o when 9 increases from 15° to 75° (see lower right panel 
of Figure [8]). On the other hand, for the case a = 0.998 M, r^k = 10 r g , and I„ oc r~ 3 , Rx/o 
increases by a factor of 10 when 9 increases from 15° to 75°. For the nearly edge on case, 
Rx/o can increase by another factor of 10 in extreme cases with small X-ray emission sizes 
and steep profiles. 

We plot the observed monochromatic flux modulated by strong lensing and geometrical 
projection as a function of inclination angle 9 in Figure [9j For the case of an extreme Kerr 
black hole case (a = 0.998 M) with steep radial profile, the strong lensing effect overwhelms 
geometrical projection, and consequently, the observed flux F v increases with 9. A comparison 
between the top and bottom panel of Figure [9] shows that difference between different spins 
is more significant for the steeper radial profile case (i.e.,n=3). 

5.2. An X-ray Ball Model 

As the second geometrical model for the X-ray emission, we assume that the X-ray 
corona is a ball around the central black hole. Such a corona is often used in the literature 
(e.g., Schnittman & Krolik 2010) to study the X-ray polarization of AGN. We furthermore 
assume that the accretion disk is opaque, and integrate over the half sphere above the 
accretion disk to obtain the observed X-ray flux. We test this X-ray ball model for three 
radii: r ball = 3r g , 5r g , and 10 r g . For the 4- velocity of the X-ray sphere, we use Eqs. (jHD 
and (II ip (Ruszkowski & Fabian 2000). We do not consider the radial motion of the X-ray 
source, and the X-ray sphere is differentially rotating around the polar axis. We show Q(9) 
for different rbaii and spins in Figure |3j Smaller rbaii gives larger Q. For small to moderate 
inclination angles, Q is dominated by u, i.e., the frame dragging, and therefore a smaller a 
gives a smaller Q, whereas for a large 9 (near the accretion disk), Q is dominated by the 
Keplerian flow Qk, and a smaller a gives a larger Q for a prograde flow, see Eq. ( TTUj) and 
Figure |2j 

We show the lensing distorted and redshifted image of the (half) X-ray ball above the 
equatorial plane in Figure [TU] for an observer at 9 = 75°, for three ball sizes, rbaii = 3r 5 , 
5r 9 , and 10 r g , and three spins, a = 0, 0.5 M, and 0.998 M. The images are taken in a way 
such that the fractional area of the image of the (half) X-ray sphere with respect to the 
area of the image would have been the same for each frame, were there no gravity. The area 
amplification is most significant for the r^aii = 3r 9 case (a factor ~ 2.3 for 9 = 75°) since it is 
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closest to the black hole, and is therefore influenced by the lensing effect most significantly, 
and least significant for r ba n = 10 r g (about 20%). For the two a/0 cases, we observe that 
the image shape (not the color) is not mirror symmetric with respect to the polar axis, this 
is most obvious in the rbaii = 3r s case, because of frame dragging. The image is stretched 
more on the right hand side (the receding side). Both the gravitational redshift and the 
Doppler red/blue shift effects are more significant for smaller rbaii (see also Figure |3]). For 
example, for the moderate spin case (a = 0.5), the redshift factor g = v jv e lies between 
0.63 and 1.26 for r ball = 10 r g and between 0.22 and 1.39 for r ball = 3r g . 

For a distant observer (r Q b s ^> 1) at an inclination angle 9, the projected image area of 
a half sphere orthogonal to the line of sight is 7rr 2 (l + cos#)/2 when the space-time is flat, 
whereas the lensed image area in Kerr space-time has to be integrated numerically. Using 
Eq. f[T9~l) the monochromatic flux amplification is found to be 

Mx =^7^r- (26) 

(the l/r n factor drops out since it is a constant for the ball model, and therefore, the results 
do not depend on the radial profile). We plot M x and the lensing magnification ratio R^/o 
in Figure [TT1 (the first row). As mentioned earlier, the amplification in the image area is the 
largest for the fbaii = 3r 9 case (Figure [TUJ). However, since the redshift effect dominates over 
the area amplification, the rbaii = 3 r g case shows the smallest amplification in the observed 
flux (Figure [TTj top left panel). As for the inclination angle dependence of flux magnification 
ratio R^jq ■, we find that the magnification ratio changes by less than 5% (15° < 9 < 75°) for 

all three spins when rbaii = 10 r g . For smaller rbaii, we find that the magnification ratio R^/o 
increases by ~20-40% for different spins when r ba ii = 5 r g and by ~80-130% for r ba n = 3 r g 
with different spins. For large inclination angles, the magnification ratio Rx/o decreases 
with 9 which is different from the X-ray disk model in the previous subsection (see Figures 
[7]and|8]), and is mainly caused by the fact that in the case of no lensing the projected area 
of the thin optical disk approaches zero when 9 approaches 7r/2 (the amplification of image 
area therefore approaches infinity) whereas that of the (half) X-ray ball approaches nr 2 /2 
(the area magnification is therefore finite). 

Since the geometry of the X-ray ball model is different from that of the optical emission 
(a thin disk), the projection ratio R^q = (l + cos0)/2cos0 is not a constant. Consequently 
both Kerr strong lensing and the projection effect contribute to the inclination angle depen- 
dence of the observed X-ray-to-optical flux ratio, and Rx/o 7^ Rx/o ( see ^1- ESJ) - We plot 
Rx/o in Figure ITT1 (the bottom right panel). Rx/o can increase by a factor of ~10 for large 
inclination angle (9 > 70°) where the major contribution comes from the geometrical pro- 
jection instead of general relativity. We also show the observed flux as a function of 9 in 
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Figure [TT] (the bottom left panel). Similar to the low X-ray disk model, Fx can increase with 
9 for small ball radii. 



5.3. High X-Ray Disk Model 

As the last example of a X-ray emission model, we consider a thin X-ray disk of radius 
10 r g which is above the black hole (parallel to the equatorial plane) with height z = 3 r g , 5r g , 
and 10 r g , respectively. This is motivated by the light bending model of Fabian & Vaughan 
(2003) where the authors tried to explain the rapid variability in the continuum observed 
for MCG6— 30— 15 using the fact that gravitational light bending by the strong field of the 
central black hole amplifies small changes in the X-ray source height (see also Miniutti et al. 
2003; Miniutti & Fabian 2004). Intuitively, a larger z means greater distance of the X-ray 
disk from the black hole, and therefore less influence from the strong gravitational field, and 
this might give smaller inclination angle dependence in Rx/o when 9 is not too large. We 
plot the angular velocity profile in Figure HJ From Figure HJ the disk with the smallest height 
(closest to the black hole) has the largest velocity (for the same p, distance to the polar axis, 
and spin parameter a). 

We show images of the lensed disk above the black hole as seen by an observer at 9 = 75° 
in Figure [T21 for three heights, and three spins. The area amplification is most significant for 
the height z = 3r g case (a factor of ~1.9 for all three spins), and least significant for the 
z = 10 factor about 1.5 for all three spins). The redshift effect (gravitational and 

Doppler) is also most significant for the z = 3 r g case (see the first column of Figure fl~2l . For 
example, for moderate spin a = 0.5, the g factor lies between 0.84 and 1.02 for z^isk — 10 r g 
case, and between 0.58 and 1.16 for z disk = 3r g case. Since points in the interior of the disk 
are closer to the black hole, they are more influenced by gravitational redshift. Since the Q(p) 
profile (Ruszkowski & Fabian 2000) is not monotonic for this model (note that p increases 
with 9 for fixed Zdisk); the points on the disk with most significant Doppler blue/redshift 
can lie in the interior of the disk for small Zdisk- For example, when z^isk = 3r g , the largest 
redshift occurs in the interior of the disk instead of on the boundary (Figure [12j the first 
column) . 

We choose no cutoff at the center of the disk, and use the radius-independent intensity 
profile (since the smallest disk height we used is 3r g , the dependence of the results on the 
radial profile should be much weaker than the low X-ray disk model in §5.11 where the case 
most sensitive to the radial profile is when a = 0.998 and risco = 1.23 r s ). This gives less 
weight to the central part of the disk which is closest to the black hole compared with the 
n = 3 (l v oc r~ 3 ) case, and will yield more conservative constraints. The projection effect 
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(i.e., the cos 9 factor) is the same for X-ray and optical emission, consequently Rx/o — -^x/o • 
We show the magnification ratio Rx/o m Figure [131 Since the black hole spin contributes 
to the space-time curvature only in higher orders, and the smallest X-ray disk height z 
we test is 3r g , we found that the results for different spins do not differ significantly. For 
15° < 6 < 75°, the magnification ratio Rx/o (flux) increases by about 85%, 55%, and 
27% for z = 3r g , 5r g , and 10 r g respectively Therefore, a thin X-ray disk of radius 10 r g 
with Ruszkowski & Fabian (2000) velocity profile is consistent with observed a ox for normal 
quasars even when it is above the central black hole with a height as low as 3 r g (Rx/o ~ 2 
corresponds to Aa ox ~ 0.1, well within the observed scatter of a ox , see next section). For 
high inclination angles (corresponding to BALs), the magnification ratio Rx/o in observed 
monochromatic flux is larger than one (say, a factor of 2, see Figure [TBI which is similar to 
the low X-ray disk model (although less significant). This implies that the absorption in the 
X-ray band might be under-estimated. 

In the so-called light-bending model (Miniutti et al. 2003; Fabian & Vaughan 2003; 
and Miniutti & Fabian 2004) the authors used the variation of the X-ray source height in 
a strong gravitational field as one interpretation for the observed large variability of the 
X-ray power-law component compared to that of the reflection component (in particular the 
fluorescent iron line emission). Choosing a ring- like source with distance p s = 2r g from the 
z-axis and assuming that the source is co-rotating with the accretion disk with the same 
orbital velocity as the underlying disk, Miniutti & Fabian (2004) found that the observed 
flux from the (direct) power-law component can increase by a factor of ~20 when the ring 
hight z increases from lr g to 20 r g , or a factor of ~5 when z increases from 3r g to 10 r g (for 
inclination angle 9 = 30°). For our disk model with r^sk — 10 r g and Ruszkowski & Fabian 
(2000) velocity profile, we found that if such a X-ray disk is observed with an inclination 
angle 6 = 30°, the observed X-ray flux can increase by a factor of ~40% when the disk 
height changes from 3r g to 10 r g (see Figure [TBI . Our results are qualitatively consistent 
with Muniutii & Fabian (2004). The quantitative difference, i.e., our models shows much 
smaller variation in observed flux of the power-law component, is not surprising, since we 
have assumed a very different geometry for the source emission (a disk of size 10 r g versus a 
ring with radius 2r g ). 
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6. Comparison With Observations 

An important parameter characterizing the flux ratio between X-ray and optical emis- 
sion is a ox (Tananbaum et al. 1979) 

a ox = 0.3838 log p^-. (27) 

2500 A 

Some authors use 3000 A for the optical emission with the constant in the front as 0.372, 
and this does not change the conclusion of this paper. There is evidence showing that the 
observed a ox in quasars are strongly anti-correlated with quasar luminosities in both the 
optical and X-ray bands; however, no significant redshift evolution was found (Strateva et 
al. 2005; Steffen et al. 2006; Lusso et al. 2010). The intrinsic scatter of a ox is about 0.2- 
0.3 at a given luminosity (Steffen et al. 2006; Lusso et al. 2010). There are several factors 
contributing to this scatter including quasar variability and intrinsic differences between 
quasars. 

In this paper, we have investigated variations in Rx/o due to lensing and projection ef- 
fects which are functions of the inclination angle 9 (see Figures EJ El and [TBI, and therefore, 
the orientation effect also contributes to the distribution of a ox . To test the significance of 
strong lensing effects on the observed X-ray-to-optical flux ratios, we separate the contribu- 
tions to a ox as, 

a ox = ajf + 0.3838 hgR x/0 , (28) 

where Rx/o includes the orientation induced contribution and a™ tr represents other contri- 
butions such as spectral variability and intrinsic differences between quasars. We assume 
the intrinsic distribution of a™ tr to be Gaussian with mean value /i = —1.5 and scatter 
a = 0.2. Assuming a uniform distribution for a™ tr > e -g-> between (—2, —1), does not change 
our conclusions. The distributions of Rx/o are calculated for the three X-ray models which 
we have considered, assuming observers are randomly distributed in all solid angles, i.e., a 
uniform distribution in the variable fi = cos# (0 < fi < 1). We use Monte-Carlo simulations 
to combine the two distributions to obtain the distributions of a ox using 500,000 realizations. 

We assume that differences between normal quasars, BALs, and obscured quasars are 
pure observational effects caused by different inclination angles. We furthermore assume 
that obscured quasars are about 20% of the total population (i.e., Type I + Type II quasars, 
Hasinger 2008) since we are focusing on the high luminosity regime, and BALs count for 40% 
of Type I quasars (Dai et al. 2008; Shankar et al. 2008). Although the 20% Type II fraction 
might be lower than some other measurements, here we use this as a fiducial number to test 
the predictions of our model. In addition, some BALs such as LoBALs and FeLoBALs may 
have Compton thick absorptions, which can be classified as Type II quasars. The fraction of 
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these objects is non-negligible, 4-7% (Dai et al. 2012), which can compensate for our lower 
assumed Type II fraction. Consequently, we assume 0° < 6> Normal < 58.7° < 9 BAL < 78.5°. 
We do not calculate the observed a ox distributions for Type II quasars because their a ox 
is hard to study due to the heavy Compton-thick absorption. This also avoids very large 
inclination angles where the thin disk approximation assumed for optical emission and the 
two X-ray disk models can be inaccurate (we have assumed zero thickness for the thin disk 
models; however, the finite thickness of the disks must be considered for nearly edge on 
cases when computing unlensed quantities such as image area). The inclination dependence 
of Rx/o is the same as the strong lensing R^q for the low and high X-ray disk models, 
because the projection factor cos 8 is canceled out considering both the optical and X-ray 
emission. For the X-ray ball model, the projection effect, -Rx/o^ = (1 + cos#)/2cos#, * s n °t 
canceled out, and both lensing and projection effects contribute to the inclination dependent 
flux ratios. 

The results for the three X-ray models are shown in Table [2] and Figures [T4T - [T6l For 
each X-ray model, we generate the distribution of the observed a ox , and computed its mean, 
standard deviation, and skewness, i.e., (a ox ), a Qox , and S, respectively, through Monte-Carlo 
simulation. The extra scatter added to a ox by the Kerr strong lensing (also by different 
geometries for the X-ray ball model) is insignificant. We have assumed ct™ tr to be Gaussian 
with cr = 0.2 whereas the largest a in the observed a ox is less than 0.24 (corresponding 
to a (lcns) < 0.13) among all models considered in Table EJ Because the scatter of a ox is 
constrained in normal quasars (non-BALs, 0° < # Normal < 58.7° in our model), the differential 
lensing amplification between optical and X-ray emission is small at these small inclination 
angles. Therefore, to match the observed scatter in a ox , 0.2-0.3, we have to assume most of 
the scatter comes from other factors not related to the orientation effect, such as spectral 
variability or intrinsic differences between quasars. The skewness of the simulated a ox is also 
small with the largest case about 0.03. This is not surprising since distribution of a ox is 
the convolution of a™ tr with zero skewness and Rx/o, where we assumed a larger scatter in 
a™ tr . If the measured intrinsic scatter of a ox in normal quasars were smaller, the skewness 
of the distribution would be a powerful tool to constrain corona geometry. Given the large 
measured scatter in non-BALs, which requires a larger a 1 ^, it is probably not realistic to 
try to discriminate corona geometries based on the extra scatter or skewness of the lensed 
distribution of a ox . 

Nevertheless, we find that the mean of a ox , (a ox ), can be shifted significantly from the 
intrinsic values. For example, for the low X-ray disk model with the extreme spin and the 
steep radial profile, (a ox ) Normal is smaller than (a ox ) mtrmsic (equal to —1.5) by ~0.25. The 
difference in a ox between normal and BAL quasars (see the last column of Table [2]) can be 
as large as ~0.2. For the low X-ray disk and the X-ray ball model, the difference in (a ox ) 
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between normal and BAL quasars is > 0.1. For the high X-ray disk model, the splitting 
of (a ox ) between BALs and non-BALs is < 0.06. Therefore, comparing large samples of 
observed, absorption corrected a ox between BALs and non-BALs, it may be possible to 
differentiate corona geometry using the orientation effects. Here, the absorptions in BALs 
should be determined from spectral properties without assuming the same SEDs for BALs 
and non-BALs. 

A number of studies show that BALs and normal quasars have the same intrinsic spectral 
energy distribution, if the observed soft X-ray weakness of BALs is caused by intrinsic 
absorption with column densities ~10 22 -10 24 cm -2 (e.g., Green et al. 2001; Gallagher et al. 
2002). For all three X-ray models considered in this paper, we found (a ox ) BAL > (ct ox ) normal 
(Table [2j). This implies that on average X-ray emission is boosted more for BALs than 
for normal quasars. Therefore, the aforementioned high column densities might still be an 
underestimate of the intrinsic X-ray absorption of BALs. For example, let us consider a BAL 
quasar at z = 2 with JVh = 2 x 10 23 cm -2 determined to correct the a ox of this BAL to be 
the same as non- BALs. If Rxfo/ R x/o al = 2 between a BAL and a normal quasar based on 
our model calculations, corresponding to Aa ox ~ 0.11 and a typical value based on Table |2J 
then an absorption of column density ~10 24 cm -2 is needed for the BAL after including 
the strong lensing boost of the X-ray emission for BALs. Even the smallest shift in (a ox ) 
between BAL and normal quasars, i.e., 0.02 (corresponds to Rx/o ~1-13, see Table [2]), can 
cause an under-estimate of the intrinsic absorption by ANn = 2 x 10 23 cm -2 . It is therefore 
important to consider the effects of the Kerr strong gravity in order to accurately measure 
the X-ray absorptions in BALs. 



7. Conclusion 

We have developed a ray-tracing code for the Kerr metric based on numerical integra- 
tions of the 1 st order Hamilton equations for massless free particles in curved space-time 
using an adaptive 5 th order Runge-Kutta algorithm. This is the first step toward developing 
a fully general relativistic 3D radiation transfer code for the Kerr space-time with arbitrary 
flows using characteristic methods and operator perturbation (Baron et al. 2009b). We de- 
rive the formalism for studying the strong lensing of the optical and X-ray emission by Kerr 
black holes in §|3J We find that the redshift effect caused by strong gravity plays an impor- 
tant role in computing the observed flux. We do not consider the cosmological redshift, but 
to include the cosmic redshift is trivial: In Eqs. ( 117]) and f fl9|) . replace the strong lensing 
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redshift factor g by g/(l + z s ), where z s is the cosmo logical redshift of the quasar o 

Using the sizes of the optical and X-ray emission regions measured by recent quasar 
microlensing observations (Kochanek 2004; Pooley et al. 2006; Morgan et al. 2008; Eigenbrod 
et al. 2008; Chartas et al. 2009; Dai et al. 2010; Blackburne et al. 2011; Blackburne et al. 
2012), we investigate the effects of the strong gravity field of the central Kerr black hole 
on the optical and X-ray continuum emission. In particular, we calculate the correction to 
the X-ray-to-optical flux ratio caused by differential lensing distortions for the X-ray and 
optical emission. We assumed the standard thin disk model (Shakura & Sunyaev 1973) for 
optical emission and find that the effect of the strong gravity is not very significant on the 
optical emission, since its emission size, ~100r 9 , is about 10 times bigger than the X-ray 
emission (Figure E]). The effect of strong gravity on optical emission is important only for 
large inclination angles. We assume a simple power law and test three simple geometries for 
the X-ray emission: an X-ray disk immediately above the accretion disk with inner cutoff 
at risco moving with the Keplerian flow, an X-ray ball surrounding the central black hole 
moving with Ruszkowshi & Fabian (2000) velocity flow, and an X-ray disk above the black 
hole with nonzero height moving with Ruszkowshi & Fabian (2000) velocity flow. We find 
that the redshift effect is more important for the X-ray emission than for the optical emission 
(see e.g., Eqs f|T7|) and f lT9|) . Figures [5] and [11]) . Unlike the optical emission, the effect of Kerr 
strong lensing on the observed X-ray flux can be important even for moderate inclination 
angles (i.e., for normal quasars). Among the three X-ray models tested in this paper, the low 
X-ray disk model shows the largest variation in the X-ray-to-optical magnification ratio Rx/o 
(as large as a factor of 10 for normal quasars, and another factor of 10 for BAL+obscured 
quasars). Scatter by a factor of 10 (corresponds to Aa ox = 0.384) is still marginally consistent 
with current observations of a ox (see Eq. (3) of Steffen et al. 2006) for normal quasars. 
The strong inclination angle dependence makes Rx/o a random variable with nontrivial 
distribution, and the observed distribution of a ox is therefore modulated by that of Rx/o- 
In particular, a ox for a BAL can be greater than that of a normal quasar by ~0. 1-0.2 (see 
Table [2]). Intrinsic absorption with column density ~10 24 cm -2 is needed in order to balance 
this boost of X-ray emission by Kerr strong lensing. Comparing with current estimates of 
the BAL absorption column density (~ 10 22 -10 24 cm -2 , e.g., Green et al. 2001; Gallagher et 
al. 2002), we find that it is important to include the Kerr strong gravity in order to correctly 
quantity the BAL absorption. 

We note that for normal quasars, since the strong lensing distortion in the optical 



4 Note that for the optical emission, the 1 + z s factor is absorbed into the definition of the scale radius, 
i.e., the r s is the emission size at frequency v — (1 + z s )v . For X-ray emission, there will be a constant 
factor (1 + z s )~ (r+2) in front of the integral in Eq. (fT9l) . 
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emission is not significant given the larger sizes of the optical emission region (of order 
10 2 r g ), the current microlensing results need only minor modifications in order to incorporate 
the effect of strong gravity. As can be seen from Figure \5[ the strong lensing effect is 
significant for optical emission only when the inclination angle 9 is high (observed nearly 
edge on). Abolmasov & Shakura (2012) studied three high amplification events from optical 
observations of SBS J1520+530, and Q 2237+0305 using straight line caustic crossing, and 
found that a model incorporating the strong lensing of the central black holes improves the 
fit; however, it requires high inclination angles (6 > 70°) not typical for normal quasars. 
This is qualitatively consistent with our results (see §1]). 

For the X-ray emission, if the size of the X-ray emitting region is of order 10 r g (or even 
smaller, given the fact that most current microlensing observations give no lower limit on 
the X-ray emission sizes) and it is centered around the central black hole, then the X-ray 
emission will be strongly lensed by the central black hole before it is microlensed by stars in 
a lens galaxy much farther away. Therefore, it might be important to add a strong lensing 
piece to the microlensing code which in principle should give more accurate constraints on 
the source sizes, in particular, the X-ray emission sizes. Furthermore, since Kerr strong 
lensing depends on important parameters such as the spin of the black hole, the inclination 
angle of the accretion disk, and the geometry of the corona, it is possible to extract these 
valuable information by combining the Kerr and microlensing models. A first attempt along 
this direction will be presented in a companion paper (Chen et al. 2012b). 

Of course, models such as a razor-thin X-ray disk above the accretion disk with an 
isotropic intensity profile may be too simplistic. More complicated and physically solid 
models should be tested in the future. But these simple models still show the importance of 
studying the effect of the strong gravity if the X-ray source is indeed very small in size and 
very close to the central black hole. Although we have used the innermost stable circular 
orbit as a natural cutoff for the optical emission and the first X-ray model, it is possible 
to extend the emission region to within rjsco with the gas following the so-called plunging 
trajectories along geodesies (Krolik 1999; Agol & Krolik 2000) and strong gravity plays 
an even more important role — the effects shown in this paper might be more significant. 
A radiation transfer code based on the ray-tracing code developed in this paper is being 
developed by the PHOENIX group and will be published in the future §. 

We thank Ronald Kantowski for helpful discussions. We also thank the anonymous 



5 PHOENIX is a general ized radiative transfer cod e that has been developed into 3D Fra m ework in a 
series of recent papers (see lHauschildt fc Baronl l2006t Ichen et al. II2OO7I baron et al. Il2009al |bl bold and 
references therein) that can handle radiative transfer in arbitrary metrics, with arbitrary velocity fields. 



-23 - 



referee for the careful review of this work. BC and XD acknowledge support for this work 
provided by the National Aeronautics and Space Administration through Chandra Award 
Number GO0-11121B, G01-12139B, G02-13132A issued by the Chandra X-ray Observatory 
Center, which is operated by the Smithsonian Astrophysical Observatory for and on behalf 
of the National Aeronautics Space Administration under contract NAS8-03060. BC and XD 
acknowledge support for program number HST-GO-11732.07-A provided by NASA through 
a grant from the Space Telescope Science Institute, which is operated by the Association 
of Universities for Research in Astronomy, Incorporated. BC and EB acknowledge NSF 
AST-0707704, and US DOE Grant DE-FG02-07ER41517 and support for program number 
HST-GO-12298.05-A provided by NASA through a grant from the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy, 
Incorporated, under NASA contract NAS5-26555. 



A. Appendix 

We give the equations of motion (six 1 st order ordinary equations) used in our the 
backward raytracing code. Inserting Eq. (j3]) into Eq. the equations of motion follow 
after straightforward but tedious algebra. For the 4 configuration-space variables (t,r,9,<p), 
we have 

dt p t 2 Mar 

dX = ~^~^K n 

dr A 

dX = 7 Pr 

d9 pe 

dX = 7 

A-a 2 sin 2 fl\ 2Mar 



i\ \ - \ • ->,j lP* 2A P^' ( Al ) 

dX \ (rAsm 9 J p z A 

and for the two momentum-space variables p r , and pg (pt and p^ are constants of the motion), 



1 8 ( 1 \ 2 d ( 2Mar \ 1 d ( A\ 2 1 d 

2d~r\^) Pt + d~r VTA] Vm ~ 2d~r Pr ~ 2d~r 1 ~ 1 P ° 

1 d (A -a 2 sin 2 s 



2dr \ p 2 Asin 2 # 



(A2) 



and 



dX 2 89 WJ Pt + 86 V P 2 A J Pm 2 89 \ P 2 ) Pr 2 89 \p 2 



1 8 (A - a 2 sin 2 fl 
289 V p 2 Asin 2 ^ 



(A3) 



where 



and 









(a 2 


i.. 




de 1 





2M\ Pa 4 -r 4 2r 2 a 2 sin 2 fl 



Ap 2 y V A P 2 

AMa 2 r(a 2 + r 2 ^ 



Ap 4 



sin e cos 0, 



— (-) = l(r-M- r 4 
dr \p 2 J P 2 \ P 2 

d ?A\ 2 2 

-^-j = -aAsintfcostf, 



d ( 1 \ ^ r 



<9r Vp 2 / P 4 

^ ^ 2 ■ n n 

de\7) = V s1 cos ' 



d (2Mar\ _dg 03 _ 2Ma [p 2 {a 2 - r 2 ) - 2r 2 A] 

dr~ V P 2 A ) ~ dV ~ pW 
d (2Mar\ dg 03 _ AMra 3 sin e cos 9 

W\ p 2 A ) ~ ~~df~ p^A ' 



d (A- a 2 sin 2 9 \ dg 33 2a 2 sin 2 9 [r A + (r - M)p 2 ]- 2r A 2 



<9r V p 2 Asin 2 # J dr p 4 A 2 sin 2 9 

d (A- a 2 sin 2 0\ dg 33 _ -2 cos #[p 2 (p 2 - 2Mr) + 2Mra 2 sin 2 9} 



de \ p 2 Asin 2 # y de p 4 Asin'^ 
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Table 1. Strong Lensing Flux Magnification for Optical and X-ray Emission. 



Model a Inclination Angle (deg) 



1° 15° 30° 45° 60° 67.5° 70° 75° 80° 82° 85° 89° 



optical emission b 

a = 0.998, r disk = 200 r g 097 097 098 099 L02 L04 L05 L09 L15 L19 L31 2A5 

a = 0.5, r disk = 200r 9 0.97 0.97 0.98 0.99 1.02 1.04 1.05 1.08 1.13 1.17 1.27 2.25 

a = 0, r disk = 200 r 9 0.97 0.97 0.98 0.99 1.02 1.04 1.05 1.07 1.12 1.16 1.26 2.17 



low X-ray disk 



a 


= 0.998, n 


= 3, r disk 


= 5r 3 


0.05 


0.06 


0.09 


0.17 


0.38 


0.63 


0.76 


1.15 


1.96 


2.57 


4.35 


22.79 


a 


= 0.998, n 


= 0, r disk 


= 5r 9 


0.15 


0.17 


0.23 


0.38 


0.71 


1.03 


1.18 


1.58 


2.27 


2.73 


3.96 


15.13 


a 


= 0.998, n 


= 3, r disk 


= 10 r g 


0.11 


0.12 


0.15 


0.24 


0.46 


0.70 


0.83 


1.20 


1.94 


2.49 


4.10 


20.62 


a 


= 0.5, n 


— 3, 7* d i sk 


= 10 r g 


0.34 


0.36 


0.44 


0.62 


0.94 


1.20 


1.30 


1.59 


2.03 


2.31 


3.04 


9.54 


a 


= 0, n : 


= 3, r disk 


= 10 r a 


0.41 


0.44 


0.52 


0.69 


0.98 


1.20 


1.29 


1.53 


1.89 


2.12 


2.74 


8.25 


a 


= 0.998, n 


= 0, r disk 


= 10 r g 


0.39 


0.41 


0.49 


0.65 


0.92 


1.14 


1.23 


1.47 


1.85 


2.10 


2.77 


8.79 


a 


= 0.5, n 


= 0, r disk 


= 10 r g 


0.42 


0.44 


0.53 


0.69 


0.97 


1.19 


1.27 


1.50 


1.85 


2;. 07 


2.66 


7.94 


a 


= 0, n : 


= 0, r diak 


= 10r 9 


0.44 


0.47 


0.55 


0.72 


0.99 


1.19 


1.28 


1.49 


1.83 


2.04 


2.61 


7.71 



X-ray ball 



i = 0.998, r ba il = 


3r 9 


0.22 


0.22 


0.25 


0.29 


0.37 


0.42 


0.44 


0.47 


0.51 


0.52 


0.55 


0.58 


a = 0.5, r ball = 


3r 9 


0.14 


0.15 


0.17 


0.20 


0.27 


0.32 


0.35 


0.40 


0.46 


0.50 


0.54 


0.62 


i = 0, r ba il = 


3r 9 


0.10 


0.10 


0.11 


0.13 


0.17 


0.20 


0.22 


0.25 


0.31 


0.34 


0.40 


0.57 


a = 0.998, r ba ii = 


5r 9 


0.45 


0.46 


0.49 


0.53 


0.59 


0.62 


0.64 


0.66 


0.69 


0.70 


0.71 


0.74 


a = 0.5, r ball = 


5r 9 


0.42 


0.43 


0.45 


0.50 


0.56 


0.60 


0.62 


0.65 


0.68 


0.69 


0.71 


0.74 


« = 0, rbaii = 


5r 9 


0.39 


0.40 


0.43 


0.48 


0.55 


0.60 


0.62 


0.65 


0.69 


0.70 


0.73 


0.76 


i = 0.998, r ball = 


10r fl 


0.72 


0.73 


0.73 


0.75 


0.78 


0.80 


0.81 


0.82 


0.84 


0.84 


0.85 


0.86 


i = 0.5, rbaii = 


10 Tg 


0.71 


0.72 


0.72 


0.74 


0.77 


0.79 


0.80 


0.82 


0.83 


0.84 


0.85 


0.86 


a = 0, r ba ii = 


10 Tg 


0.71 


0.71 


0.72 


0.74 


0.77 


0.79 


0.80 


0.82 


0.83 


0.84 


0.85 


0.87 


high X-ray disk 


i = 0.998, z dis k = 


3r 9 


0.50 


0.52 


0.56 


0.64 


0.78 


0.90 


0.96 


1.11 


1.40 


1.61 


2.22 


8.67 


i = 0.5, z disk = 


: 3r 9 


0.49 


0.50 


0.55 


0.63 


0.77 


0.89 


0.94 


1.10 


1.40 


1.61 


2.24 


8.84 


1 = 0., Zdisk = 


: 3r 9 


0.47 


0.49 


0.53 


0.61 


0.76 


0.88 


0.94 


1.10 


1.40 


1.62 


2.27 


9.09 


i = 0.998, z dis k = 


5r 9 


0.58 


0.59 


0.62 


0.68 


0.79 


0.89 


0.94 


1.07 


1.34 


1.54 


2.13 


8.35 


i = 0.5, z disk = 


: 5r 9 


0.57 


0.58 


0.61 


0.67 


0.78 


0.88 


0.93 


1.07 


1.34 


1.54 


2.13 


8.44 


1 = 0., Zdisk = 


: 5r 9 


0.57 


0.58 


0.61 


0.67 


0.77 


0.87 


0.92 


1.06 


1.34 


1.54 


2.15 


8.55 


i = 0.998, z disk = 


= 10 r a 


0.71 


0.71 


0.73 


0.76 


0.84 


0.91 


0.95 


1.06 


1.27 


1.44 


1.93 


7.14 


i = 0.5, z disk = 


-- 10 r 9 


0.70 


0.71 


0.73 


0.76 


0.83 


0.91 


0.94 


1.05 


1.27 


1.43 


1.92 


7.15 


1 = 0., Zdisk = 


10 r 9 


0.70 


0.71 


0.72 


0.76 


0.83 


0.90 


0.94 


1.05 


1.27 


1.43 


1.93 


7.17 



a To obtain the inclination-dependent flux profile (sec Figuresl5"ll9"l llllll3jl . the geometrical projection effect should be included, 
i.e., a cos 8 factor for optical emission and the two X-ray disk models, or a (1 + cos(9)/2 factor for the X-ray ball model. 

b Mo = 0.2 + 0.8M<2oo r g , where M<200 r g is the strong lensing magnification of the optical flux within r d i sk = 200 r g . 
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Table 2. a ox Distributions After Kerr Strong Lensing. 



Model a 


Normal QSO b 




BALQSO 


b 


Normal+BAL QSO b 


A(a ox ) c 




Kx> c 


a 


s 


("ox) 


(7 


S 


(«ox) 


17 


5 




low X-ray disk 


a = 0.998, n = 3 


-1.76 


0.211 


0.003 


-1.55 


0.208 


0.004 


-1.67 


0.234 


0.03 


0.21 


a = 0.5, n = 3 


-1.60 


0.205 


0.001 


-1.47 


0.202 


0.000 


-1.55 


0.214 


0.002 


0.13 


a = 0, n = 3 


-1.58 


0.204 


0.001 


-1.47 


0.202 


0.000 


-1.53 


0.210 


0.001 


0.11 


a = 0.998, n = 


-1.59 


0.204 


0.001 


-1.48 


0.202 


0.001 


-1.54 


0.210 


0.002 


0.11 


a = 0.5, n = 


-1.57 


0.203 


0.001 


-1.47 


0.201 


0.000 


-1.53 


0.209 


0.001 


0.10 


a = 0, n = 


-1.57 


0.203 


0.001 


-1.47 


0.201 


0.000 


-1.53 


0.208 


0.001 


0.10 


X-ray ball d 


a = 0.998, rwi = 3r 9 


-1.68 


0.204 


0.000 


-1.54 


0.205 


0.003 


-1.63 


0.216 


0.014 


0.14 


a = 0.5, r ball = 3r 9 


-1.75 


0.205 


0.000 


-1.58 


0.207 


0.006 


-1.68 


0.222 


0.029 


0.17 


a = 0, r ball =3r, 


-1.82 


0.204 


0.001 


-1.66 


0.208 


0.007 


-1.75 


0.221 


0.035 


0.16 


a = 0.998, rball = 5r 9 


-1.58 


0.202 


0.001 


-1.47 


0.203 


0.002 


-1.54 


0.209 


0.011 


0.11 


a = 0.5, rbaii = 5r 9 


-1.59 


0.202 


0.001 


-1.48 


0.204 


0.002 


-1.55 


0.210 


0.013 


0.11 


a = 0, r ball = 5 r 9 


-1.60 


0.203 


0.001 


-1.48 


0.204 


0.002 


-1.55 


0.212 


0.015 


0.12 


high X-ray disk 


a = 0.998, 2 = 3r 9 


-1.58 


0.201 


0.000 


-1.52 


0.201 


0.000 


-1.55 


0.203 


0.001 


0.06 


a = 0.998, 2 = 5r 9 


-1.57 


0.201 


0.000 


-1.52 


0.201 


0.000 


-1.55 


0.202 


0.001 


0.05 


a = 0.998, 2 = 10 r 9 


-1.54 


0.200 


0.000 


-1.52 


0.200 


0.000 


-1.53 


0.201 


0.000 


0.02 



a For the X-ray ball model, we test only for case rball = 3r 9 , and 5r 9 considering the relative flatness of i?x/o(^) f° r 
r ball = 10r 9 (Figurc llll l. For the high X-ray disk model, we consider only the a = 0.998 case because of the similarity of 
the Rx/oW curves between different spins (Figure ll3H . 

b We assume that obscured (Type II) quasars are about 20% of quasar population and BALs count for 40% of Type I 
quasars. The intrinsic distribution of o ox is assumed to be Gaussian with fi = —1.5 and a = 0.2. (a x), c, and S are 
respectively the mean, standard deviation, and skewness (fi^/a 3 ) of a ox distribution. 

c A(a ox ) = (a x) BAL — (ctox ) Normal , the difference in mean of a ox between BALQSOs and normal quasars. 

d There exists extra scatter of a ox (besides strong gravity) due to projection effect of different geometries of X-ray (half 
sphere) and optical emission (thin disk). We have included this effect into our Monte-Carlo simulation. 
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Fig. 1. — Schematic plot for backward raytracing. 9 is the inclination angle, rjsco is 
the innermost stable circular orbit, uj is the dragging of the inertial frame. The dashed 
parallelogram is the intersection of the pencil of light beam (shooting backward from the 
observer) with the equatorial plane whose projection along the line of sight is the image 
window. The observer is sufficiently far away from the black hole, such that the space-time 
is nearly flat, and is close enough such that the observer has the same cosmo logical redshift. 
The accretion disk size is greatly exaggerated with respect to the observer distance r b s - A 
photon is emitted with frequency u e (measured in the local rest frame) from the disk and 
observed with gravity and Doppler red/blue shifted frequency v Q along direction n measured 
by the so-called zero-angular-momentum-observer. The light path is bent by the strong 
gravity produced by the central supermassive black hole. 
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Fig. 2. — Angular velocity Q for the Keplerian flow (see Eq. [TO]) as a function of the Kerr 
radial coordinate r for the low X-ray disk model with radius raisk = 10 r g and inner cutoff at 
?"isco = 6r 9 , 4.23 r g , and 1.24r g for spin a = 0, 0.5 M, and 0.998 M, respectively. 
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Fig. 3. — Angular velocity Q as a function of the inclination angle 6 for the X-ray ball model 
with r ba n = 3r g , 5r g , and 10 r g , and spin a = 0, 0.5 M, and 0.998 M. 
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Fig. 4. — Angular velocity Q as a function of p, distance to the polar axis, for the high 
X-ray disk model with disk height Zdisk = 3r 9 , 5r g , and 10 r g , and spin a = 0, 0.5 M, and 
0.998 M. 
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Fig. 5. — Inclination angle dependence of strong lensing of optical emission. In the top panel, 
we plot the lensing magnification Mq of image area (blue) and monochromatic flux (red). 
In the bottom panel we plot the lensed flux including the projection effect (i.e., the cos# 
factor) (normalized to 1 when observed face on). The effect of Kerr lensing is insignificant 
except for large inclination angles. We assume the optical emission is Planckian with half 
light radius r ha i f = 100 r g (scale radius r s = 41 r g ), and disk size r disk = 200 r r 
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Fig. 6. — Lensed and redshifted images for a (low) X-ray disk of radius r^k = 20 r g moving 
with the Keplerian flow. The top, middle, and bottom rows show the results for spins 
a = 0.998 M, 0.5 M, and (corresponding to inner cutoff at risco = l-24r ff , 4.23 r g , and 
Qr g ). The first, second, and last columns are for inclination angles 9 = 15°, 45°, and 75°, 
respectively. The observer is at r Q b s = 10 6 r g . The colorbar shows the redshift factor g = v j v c 
for inclination angles 9 = 15° (the first row), 45° (the second row) and 75° (the last row), 
respectively. A disk observed at a large inclination angle spans a larger redshift interval. 
The source on the left hand side is moving toward the observer (Doppler blueshifted), and 
the source on the right hand side is receding from the observer (Doppler redshifted). The 
image of an optical emitting disk is similar (but larger) since they are both in the equatorial 
plane and move with the Keplerian flow. 
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Fig. 7. — Inclination angle dependence of the strong lensing magnification of the image area 
of an X-ray disk. We assume that X-rays are emitted from a thin disk immediately above the 
disk with r disk = 10 r g and r inner = ri SC o for a = 0, 0.5 M, and 0.998 M. For the a = 0.998M 
case, we also considered rdisk = 5r 9 (red dashed line). 
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Fig. 8. — The inclination angle dependence of the strong lensing magnification of the X-ray 
flux for the low X-ray disk model. The first and second rows plot the magnification Mx and 
the X-ray-to-optical magnification ratio Rx/o, respectively. The projection effect (the cos 6* 
factor) cancels out for the X-ray disk model, and consequently the inclination-dependence 
of Rx/o is purely general relativistic. Rx/o depends on the inclination angle 9, and the 
dependence is stronger for smaller X-ray disks, larger spins, or steeper radial profiles. For 
the case a = 0.998 M, r<iisk = 10 r g , and l v oc r -3 , Rx/o increases by a factor of 10 when 9 
increases from 15° to 75°. For the nearly edge on case, Rx/o can increase by another factor 
of 10 in extreme cases with small X-ray emission sizes and steep profiles (see the bottom left 
panel). 
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Fig. 9. — The inclination angle dependence of the observed X-ray flux (normalized to 1 when 
observed face on) for the low X-ray disk model. For extreme Kerr black hole with steep radial 
profile (i.e., a = 0.998 M and n = 3, see the red curves in the top panel), the strong lensing 
magnification overwhelms the demagnification caused by geometrical projection for large 
inclination angles, and consequently the observed flux increases with 9. 
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Fig. 10. — Lensed and redshifted images for the X-ray ball model. The top, middle, and 
bottom rows show results for spins a = 0.998 M, 0.5 M, and 0. The first, second, and last 
columns are for r ball = 3r g , 5r g , and 10 r g , respectively. The observer is at r obs = 10 6 r g 
with an inclination angle 75°. The area amplification is more significant for smaller ball radii 
(a factor ~2.3 for rbaii = 3r 9 case). The source on the left hand side is moving toward the 
observer, and the source on the right hand side is receding from the observer. 
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Fig. 11. — The inclination angle dependence of the observed X-ray-to-optical flux ratio 
assuming X-rays emitted from a (half) ball around the black hole with r ba n = 3r g , 5r g , 
and 10 r g , respectively. We plot the strong lensing magnification of the X-ray flux (the top 
left panel), X-ray-to-optical magnification ratio -Rx/cf (^ ne right panel), observed X-ray 
flux (normalized to 1 when 9 = 0°, the bottom left panel), and the total correction ratio 
Rx/o = ^x/o^x/o^ (^ ne bottom right panel). 
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Fig. 12. — Lensed and redshifted images for the high X-ray disk model. The disk size 
is r disk = 10 r g . The top, middle, and bottom rows show results for spins a = 0.998 M, 
0.5 M, and 0. The first, second, and last columns are for height z disk = 3r 9 , 5r g , and 10 r g , 
respectively. The observer is at r G b s = 10 6 r s with an inclination angle 75°. The source on 
the left hand side is moving toward the observer, and the source on the right hand side is 
receding from the observer. 
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Fig. 13. — Inclination angle dependence of the X-ray-to-optical flux ratio for the high X-ray 
disk model. We assume that X-rays are emitted from a disk with r^sk — 10 r g , above the black 
hole, with height z = 3r g , 5r g , and 10 r g , respectively. We plot the lensing magnification of 
the X-ray flux Mx (the top left panel), the observed flux (normalized to 1 when observed 
face on, the top right panel), and the X-ray-to-optical magnification ratio Rx/o — Rx/o 
(the bottom panel) as function of inclination angle 9. 
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Fig. 14. — Histograms of a ox for normal (blue) and BAL quasars (green) including the Kerr 
strong lensing effects for the low X-ray disk model. The intrinsic scatter of a ox is assumed to 
be Gaussian with /i = —1.5, a = 0.2 (dashed red line). (a ox ) Normal is smaller than (a ox ) BAL 
by -0.1-0.2. See also Tabled 
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Fig. 15. — Histograms of a ox for normal (blue) and BAL quasars (green) including the Kerr 
strong lensing and projection effects for the X-ray ball model, where the extra projection 
effect is caused by the different geometries of X-ray and optical emission. The intrinsic 
scatter of a ox is assumed to be Gaussian with fi = —1.5, a = 0.2 (dashed red line). (a ox ) for 
normal quasars is smaller than that of BALs by ~0.15 and ~0.11 for rbaii = 3r 5 , and 5r g , 
respectively. The contribution of the projection effects to A(a ox ) is ~0.08, independent of 
the size of the X-ray ball. See also Table EJ 
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Fig. 16. — Histograms of a ox for normal (blue) and BAL quasars (green) including the Kerr 
strong lensing effects for the high X-ray disk model. The intrinsic scatter a ox is assumed to 
be Gaussian with \i = -1.5, a = 0.2 (dashed red line). We find (a ox ) BAL - (a ox ) Normal < 0.06. 
See also Table [2j 



